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ABSTRACT 

Precise exoplanet characterization requires precise classification of exoplanet host stars. The masses 
of host stars are commonly estimated by comparing their spectra to those predicted by stellar evolution 
models. However, spectroscopically determined properties are difficult to measure accurately for 
stars that are substantially different from the Sun, such as M-dwarfs and evolved stars. Here, we 
propose a new method to dynamically measure the masses of transiting planets near mean-motion 
resonances and their host stars by combining observations of transit timing variations with radial 
velocity measurements. We derive expressions to analytically determine the mass of each member 
of the system and demonstrate the technique on the Kepler-18 system. We compare these analytic 
results to numerical simulations and find the two are consistent. We identify eight systems for which 
our technique could be applied if follow-up radial velocity measurements are collected. We conclude 
this analysis would be optimal for systems discovered by next generation missions similar to TESS or 
PLATO, which will target bright stars that are amenable to efficient RV follow-up. 

Subject headings: Methods: analytical; Planets and satellites: fundamental parameters; Stars: funda- 
mental parameters 



1. INTRODUCTION 

With modern radial velocity techniques and the phe- 
nomenal success of space-based transit surveys, exo- 
planetary science has moved from a "stamp-collecting" 
era of finding individual systems to an era where hun- 
dreds of planetary sy stems are discovered simultaneously 
(|Borucki et al.l l20Tlf ) . Despite these successes, accurate 
characterization of planets still presents a challenge. In 
general, uncertainties in the radii and masses of planets 
are dominated by unce rtainties in the radii a nd masses 
of their host stars (e.g. iMuirhead et al.1 12012). Difficul- 
ties in characterizing the physical properties of planets 
are particularly acute for systems discovered by the Ke- 
pler space telescope. For many systems, the ratio be- 
tween the radius of the planet and the ra dius of its host 
star is known to within 1 part in 1000 ( ]Batalha et al.l 
2012). Yet the stellar radii are often not known to bet- 
ter than ten percent, meaning much of the precision 
of Kepler is lost when estimating planeta ry properties 
(jLissauer et al.ll2012l : Dohnson et alj|2012bl ). 

In general, measuring the masses of exoplanet host 
stars is a model-dependent procedure. For nearby stars 
with trigonometric parallaxes, one compares the lumi- 
nosity, effective temperature, a nd metallicity of a star 
to stellar evolution m odel grids (jValenti fc Fischenl2005l : 
Uohnson et al.|[2012al ) . For stars without measured paral- 
laxes, the stellar density can be measured from the tran- 
sit light curve and used in place of the luminosity. How- 
ever, this relies on either the assumption that the planet's 
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orbit is circular — a poor assumption for periods larger 
than 10 days — or an RV orbita l solution (jSozzetti et al.l 
l2007t iDawson fc Johnsonl l2012f ) . The atmospheres and 
interior structures of stars are also poorly understood 
for stars that differ substantially from the Sun, compli- 
cating their analyses further. Thus, model-independent 
methods of measuring stellar masses are extremely valu- 
able; 

lAgol et all (|2005f ) suggest that in a system with tran- 
siting planets, a precise measurement of the transit du- 
ration, which depends on stellar density, coupled with 
radial velocity information and precise measurements of 
the scatter in transit times can provide a unique mea- 
surement of the stellar mass. Unfortunately, this strat- 
egy requires precise knowledge of the inclination of the 
system, which from a transit lig ht curve is deg enerate 
with limb-darkening coefficients (|Jha et al.ll2000D . espe- 
cially for low signal-to-noise transit detections. 

The method described by Agol et al. also breaks down 
for resonant systems, as it assumes the relative positions 
of the planets change from transit to transit. Moreover, 
outside of resonance, transit timing effects arc small for 
all but the largest planets, so this method is not ideal for 
studying rocky planets. This strategy is successful when 
the perturbing obj ect is massive, as is the case in cir- 
cumbinary planets (|Dovle et al.ll20lH iWelsh et al.1l2012f ) 
but is less promising for studying solar-type systems. It 
has also been suggested that in a system containing a 
transiting planet and an exomoon detected through tran- 
sit timing and duration variations, the stellar mass and 
radiu s can be determ ined directly through dynamical ef- 
fects ([Kipping! I2010f ). While this technique holds future 
promise, the exomoons required to test this procedure 
have not yet been detected. 

Recently, transit timing variations caused by mu- 
tual gravitational interactions of bodies in multiple- 
planet systems have been detected (jHolman et al.ll2~010l ; 



2 



Montct and Johnson 



iFord et al]|2012al) . These deviations from a linear transit 
ephemeris allow for an estimate of the ratio of the mass 
of the perturbing planet to the mass of its star. In cases 
where multiple planets transit, the ratio of the masses of 
each planet to the ma s s of the host star can be estimated 
(|Fabrvckv et alJl2012t iSteffen et al.ll2012aft . 

In this paper, we propose a method to directly mea- 
sure stellar and planetary masses for multi-transiting sys- 
tems by combining an analysis of the transit timing sig- 
nal caused by planet-planet interactions with Dopplcr 
radial veloc ity measurements . Unlike the technique de- 
veloped bv lAgol et all (|200l . our method requires the 
observed transiting planets to lie near a mean-motion 
resonance, where transit timing effects are strongest. In 
§2, we explain how transit timing variations can be com- 
bined with radial velocity information to estimate stellar 
and planetary masses. In §3, we apply our process to 
the well-studied Kepler-18 planetary system, and com- 
pare the result to both numerical integrations of the sys- 
tem and published stellar evolution models. We find the 
scheme is viable, but at present the radial velocity data 
is insufficient to place meaningful constraints on stellar 
parameters. In §4, we discuss uncertainties and limita- 
tions to our method, as well as its applications to systems 
discovered by Kepler and its eventual successors. 

2. UNIQUE MASSES AND ERRORS 

2.1. Mass Determination from TTVs 

Consider a system of two coplanar planets orbiting 
near (but not exactly at) a first-order mean motion res- 
onance. The planets have periods P and P' (here and 
throughout, the unprimed quantity refers to the inner 
planet and the primed quantity to the outer planet) and 
orbit their star such that the inner planet completes ap- 
proximately j orbits in the time the outer planet com- 
pletes j — 1. A nearly edge-on observer will detect both 
planets transiting their host star. Because of the near- 
commensurability of their periods, the inner planet will 
pass its companion at nearly the same location each or- 
bit, driving small gravitational interactions which add 
coherently, inducing a small forced eccentricity on each 
object. The two planets will therefore not transit their 
star in an exactly periodic fashion. Instead, a small, sinu- 
soidal departure from periodicity, termed a transit timing 
varia tion (TTV), will be observed (e.g. iNesvornv et al.l 
12012ft . TTVs have be en used to detect the presence of 
nontransiting planets (Ballard e t al.l[20TTI : iDawson et al.1 
l2012f ) and to fully characterize systems when multi- 
ple p lanets transit (jHolman et al.1 120101 : iLissauer et al.1 
120111) . The period of the TTV signal is related to the 
periods of the planets such that 
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In most cases, the superb photometry provided by the 
Kepler mission allows this quantity to be precisely esti- 
mated. 

An analyti c form for the amplitu de of the TTV signal 
is derived bv iLithwick et al.l (|2012L hereafter LI 2). The 
amplitude of the signal depends strongly on the free ec- 
centricity of the system. Here, free eccentricity refers to 
the component of the eccentricity caused by the initial 
dynamical conditions of the system, not the component 



driven by resonant interactions. Without observing a 
secondary transit, for small planets the free eccentricity 
is difficult to constrain precisely via photometry. How- 
ever, many TTV signals have been detected in systems in 
which the planets have orbital periods of days to weeks. 

In this case, the estimated ages of the planets are larger 
than the expected tidal circularization timescale at their 
present locations, so their orbits can be expected to have 
negligible free eccentricity. This can be verified by ana- 
lyzing the phase of the TTV signal. If the zero points of 
the TTV signal occur when the longitude of conjunction 
is parallel to the line of sight, LI 2 suggest the free eccen- 
tricity can be neglected. In this case, the amplitudes of 
the TTV signals, V and V, are 
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where m and M are the planet and stellar mass, A is 
the fractional distance from commensurability, typically 
of order 0.01, and / and g the appropriate coefficient 
of the disturbing function, which characterizes the in- 
teractions between the planets. These sums of Laplace 
coefficients can be calc ulated by using the inform ation 
found in Appendix B of IMurrav fc Dermottl (|1999D . Ad- 
ditionally, the values of / and g for common resonances 
are listed in L12. To first order, these coefficients are of 
order unity and depend only weakly on A. For systems 
with TTVs, Equations [2] and [3] enable a unique determi- 
nation of the planet-star mass ratio, but normally one 
must rely on stellar models to separate stellar and plan- 
etary properties. However, if radial velocity measure- 
ments are available, the amplitude of the Doppler signal 
can be used in conjunction with the TTV information to 
estimate the masses of the planets and the star. 

2.2. Including Radial Velocities 
The scmiamplitudc of a radial velocity Doppler signal 

/ 2irG \ to sin i 
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with i and e the inclination and eccentricity, respectively 
(|Paddockl Il913f ). Despite the lack of free eccentricity, 
we may expect a small forced eccentricity as a result 
of 3-body interactions. The magnitude of this forced 
eccentricity is < 0.05, so neglecting it will induce an error 
of < 0.1% in our scmiamplitudc calculation. An error 
of the same magnitude but in the opposite direction is 
induced by assuming i = 90°, since a strong constraint 
on the inclination is provided by our requirement of a 
transit. 

Thus, neglecting eccentricity and assuming an edge- 
on orbit, in the limit where m <^ M the radial velocity 
semiamplitude can be approximated as 
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and the radial velocity can be modeled as 
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with t c the time of transit center. Again, a degener- 
acy exists between the planet and stellar mass, so stel- 
lar models must be invoked. However, the degeneracy 
is different from the one recovered from transit timing 
variations, so these two expressions taken together can 
be used to solve for the planet and stellar masses indi- 
vidually. This allows for two independent measurements 
of the mass of the star and one unique measurement of 
the mass of each planet. The mass of the star is 
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and the mass of each planet is 
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Simply put, for a given system of two coplanar planets, 
the constants M 1 / 3 V'K~ 1 and AP^VK'" 1 depend only 
on the system architecture. Thus, by precisely measuring 
the RV scmiamplitude and magnitude of the TTV sig- 
nal, the stellar mass can be directly estimated. Because 
of Kepler's exceptional photometry, the periods of each 
planet and terms derived from these (such as A) are well 
known. Thus we expect the errors in the mass estimates 
to be dominated by the errors in V and K, and neglect 
the errors caused by other terms. In this case, 
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where we expect the covariant terms to be zero since 
K and V arc independently measured quantities. Here, 
the fractional uncertainties depend quite sensitively on 
the ability to measure K and V. Typically for systems 
of multi-transiting planets, only one of these quantities 
is well measured. Therefore, we would expect to only 
weakly constrain the stellar masses at present; with more 
observations the constraints will tighten considerably. In 
the limit where m 3> m 1 , K and V are much larger than 
their counterparts and can be more easily constrained. 
Thus when one planet is substantially more massive than 
its companion, one stellar mass measurement will be con- 
siderably more precise than the other and the mass of 
the more massive planet will be better constrained than 
the less massive planet. In the case where m ~ m', both 
measurements are expected to have similar uncertainties. 



3. EXAMPLE 

Kcplcr-18 (KOI 137, KIC 8644288) is a planetary 
system containing three nearly coplanar planets with 
3.5, 7.6, and 14.9 da y periods orbiting a 0.97A/© star 
(|Cochran et al.l [20111 henceforth Cll). These planets 
(137.03, 137.01, and 137.02, or Kepler-18 b, c, and d, 
respectively) were confirmed by a combination of tran- 
sit timing and radial velocity measurements. The star 
has been observed using the Kepler short cadence mode 
nearly continuously for two years, allowing for precise 
measurements of transit times over dozens of transits. 
Moreover, 18 radial velocity measurements of this K p = 
13.5 star, where K p is the apparent magnitude in the 
Kepler bandpass, have been collected over the past three 
years by the California Planet Search team with the Keck 
1 High Resolution Echclle Spectrometer (HIRES). Thus, 
enough data exist to attempt to determine the mass of 
each member of this system dynamically. 

We first fit a limb-darkened light curve to a series of 
phase-folded transits to estimate the observable transit 
parameters, such as the impact parameter and limb- 
darkening coefficients , following the OCCU LTQUAD 
routine developed by iMandel fc Agoll (|2002f) . Because 
of the high signal to noise ratio of these observations 
and the one-minute integration times, transit parame- 
ters can be easily measured from individual transits: we 
find no significant difference in these parameters or their 
uncertainties when fitting one individual transit instead 
of fitting a phase-folded transit. 

Once the shape of the light curve is modeled, we fit 
a curve of this shape to each individual transit, allow- 
ing only the time of transit center to vary. Each indi- 
vidual transit light curve consists of over 200 in-transit 
data points, allowing for measurements of the transit cen- 
ter time to sub-minute precision. We remove from our 
dataset transits that occur simultaneously with the tran- 
sit of another planet. As expected, the transits follow a 
sinusoidal deviation from a linear ephemeris; these de- 
viations, shown in Table 21 appear to be anti-correlated 
between the two planets. 

For our method to provide meaningful mass estimates, 
the primordial (free) eccentricity of the system must be 
damped on a timescale shorter than the age of the sys- 
tem. As explained in L12, if the zero points of the transit 
timing variations occur at the times at which the longi- 
tude of conjunction of the planets is equal to or 180 
degrees, then the system is likely to have negligible free 
eccentricity. We check the phase of these transit timing 
variations by fitting each TTV curve independently to a 
sinusoid. We determine parameters of this sinusoid and 
their uncertainties by minimizing the x 2 statistic through 
a Levenberg-Marquardt algorithm. Additionally, we al- 
low for a vertical offset to the sine function in the fit 
(indicative of a miscalculated time of transit center t c ), 
and also a linear trend (indicative of a miscalculated or- 
bital period). The results of this minimization can be 
found in Table [TJ Both planets are consistent with hav- 
ing zero free eccentricity and anticorrelated TTV signals. 
From these parameters, we measure a fractional distance 
from commensurability A = —2.776 x 10~ 2 and find the 
coefficients of the disturbing function to be f = —1.251 
and g = 0.5308. The amplitude of the TTV signals for 
Kepler-18 c and d can be measured to within 3.3 and 6.8 
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percent, respectively. 

These results can then be combined with radial ve- 
locity measurements in order to uniquely constrain the 
stellar and planetary masses. Cll used 14 radial velocity 
measurements to confirm this system; we analyze these 
data plus four additional observations collected between 
1 July 2012 and 1 August 2012, all of which are provided 
in Table H 

The large uncertainties in each individual observation, 
coupled with the small number of observations relative to 
the number of observed transits, suggest our mass uncer- 
tainties will be dominated by uncertainties in the radial 
velocity semiamplitude. In fact, many different solutions 
fit the RV data satisfactorily. As an example, Cll fit 
a larger RV semiamplitude for planet d than c, despite 
the fact that they find planet c to be both more mas- 
sive and nearer the star than planet d. The analysis is 
complicated by the existence of the much smaller planet 
b orbiting inside the other two planets. In this case, we 
invoke one additional piece of information. Equations [2] 
and [3] can be combined to solve for the mass ratio of the 
resonant planets, 
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(14) 



which in this case implies m'/m = 1.22 ± 0.09, where 
m! refers to planet c and m to planet d. This can be 
applied as an additional constraint in the radial veloc- 
ity fit. In the case where <trv ^ &ttv, an equivalent 
mass ratio constraint, derived from the radial velocity 
semiamplitude ratio, can be applied to the TTV fit. 

With this additional constraint, we model the RVs as 
the sum of three sinusoids of the form of Equation^ with 
three free parameters: the semiamplitude of one of the 
resonant planets (c or d; here, we fit c), the semiampli- 
tude of the innermost planet b, and an offset term, 7. We 
find the best fitting parameters to be K c = 6.89 ± 1.40 m 
s -1 , K b = 4.18±2.14 m s" 1 , and 7 = 1.30 ±1.45 m s" 1 . 
From the mass ratio above, this implies a semiamplitude 
for planet d of = 4.52 ± 0.97 m s _1 . We now have 
enough information to estimate the stellar and planetary 
masses; these results are shown in Table [3] 

When both the RV and TTV amplitudes are measured 
without invoking the extra constraint of Equation 1 141 
two independent measurements of the stellar mass can 
be calculated, one through K and V, and one through 
K' and V. However, since our value for K' is found by 
assuming a value for K, we only calculate one indepen- 
dent measure of the stellar mass. We find a stellar mass 
of O.83±O.51M , consistent with that found by Cll. We 
find the masses of planets c and d to be 18.6 ± 11.6 and 
15.4 ± 9.5M®, respectively. 

It is somewhat disappointing that the uncertainties 
in the stellar mass are so large in this example, but 
this should be considered a shortcoming in the avail- 
able data, not in the potential of our technique. Because 
most Kepler Objects of Interest (KOIs) are considerably 
fainter than typical stars probed by radial velocity sur- 
veys, follow-up radial velocity measurements are often 
carried out only to a level necessary to confirm the plan- 
etary nature of a transiting system, if any observations 
are collected at all. Thus for most systems that exhibit 
transit timing variations, radial velocity measurements 



alone are rarely precise to within even 20%. Better con- 
straints on K arc regularly achieved for stars targeted in 
radial velocity surveys, and with more follow-up obser- 
vations these mass estimates will be greatly improved. 
This is discussed more fully in 2) 

We can confirm the validity of our method by compar- 
ing our analytic result to results obtained through nu- 
merical integrations of this system. To accomplish this 
tas k, we make use of th e Systemic Console developed 
bv IMeschiari et all ((20091 ). This program is designed to 
simultaneously fit Doppler and transit timing measure- 
ments. The Console contains several built-in integrators, 
including an eighth-order Runge-Kutta scheme employed 
in this work. The Console is not designed to enable the 
user to solve for the stellar mass as a free parameter. 
We circumvent this problem by first assuming a stellar 
mass. We fix the period and mean anomaly at BJD = 
2455128.0 so that they are consistent with values found 
in Table 7 of Cll; we then allow the planet masses, incli- 
nations, and eccentricities to vary and minimize the x 2 
of the system. Once x 2 is calculated, we vary the stel- 
lar mass slightly and repeat this procedure. With this 
technique, we can map the likelihood space in both M 
and m. As shown in Table [31 both the best fitting pa- 
rameters and their uncertainties arc consistent with the 
analytic result, suggesting that our method is viable and 
that dynamical techniques can be used in conjunction 
with our analytic result to further constrain the stellar 
and planetary parameters. 

In all cases, our uncertainties are dominated by our 
20% errors in the radial velocity semiamplitudes. The 
uncertainty in the radial velocity semiamplitude will de- 
crease considerably with more radial velocity observa- 
tions. We prove this claim by simulating observations 
placed randomly between the months of June and Octo- 
ber, when the Kepler field is visible at night. We first 
find the true radial velocity of the system at that time, 
assuming K c = 7.0 m s _1 . An statistical uncertainty a 
is randomly drawn from the observed errors in previous 
HIRES measurements, and a Gaussian random number 
is drawn from a distribution Af(0, a). The radial velocity 
measurement is shifted by an amount equal to this ran- 
dom number, and the statistical uncertainty is recorded 
as a. Finally, to simulate the effects of radial velocity 
"jitter" caused by stellar pulsations, a random number is 
drawn from 7V(0, 3 m s _1 ); this value is also added to the 
radial velocity measurement. The observations are fitted 
to a combination of sinusoids as described above, and the 
stellar and planetary masses are estimated. The frac- 
tional error in the semiamplitudes for the largest planet 
as a function of the number of observations is shown in 
Fig. [1] (solid line). 

We find that, with 30 more radial velocity observa- 
tions, the uncertainty in our calculation drops by nearly 
a factor of two, from the current 61 percent to 33 per- 
cent. To provide substantially better than 33 percent 
uncertainties without obtaining 50 radial velocity mea- 
surements, we can target a less massive star. As stated 
in this method will be optimal for stars for which 
evolutionary models are less able to constrain stellar pa- 
rameters precisely, such as F-type stars, subgiants, and 
M-dwarfs. Since M-dwarfs are less massive than their 
G-type counterparts, a given mass planet around an M- 
dwarf produces a comparatively larger RV signal. Since 
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TABLE 1 

TTV Fitting Results for Kepler-18 



Planet 


Phase (deg) 


Amplitude (min) 


Period (day) 


T c (BJD - 2454900.0) 


TTV Period (day) 


18-c 


184.5 ±4.1 


5.54 ±0.18 


7.6415716(5) 


68.4071(2) 


265.1 ±2.5 


18-d 


3.2 ± 8.8 


4.46 ±0.30 


14.858941(1) 


61.1531(1) 


265.9 ± 5.3 



TABLE 2 
Keck HIRES Relative Radial 
Velocity Measurements of 



Kepler-18. 


BJD-2440000 


RV 




17 




ms — 1 


m s 1 


15076.009 


7.750 


2 


539 


15076.927 


6.950 


2 


487 


15081.024 


8.617 


1 


214 


15082.007 


-1.007 


2 


381 


15084.984 


-7.320 


2 


977 


15318.066 


3.388 


2 


625 


15322.029 


-10.093 


2 


303 


15373.004 


12.189 


2 


150 


15403.019 


24.983 


2 


915 


15405.909 


-11.692 


2 


350 


15406.881 


0.340 


2 


195 


15413.011 


-10.788 


2 


498 


15432.970 


0.205 


2 


233 


15436.782 


-6.675 


2 


256 


16109.905 


-14.919 


2 


261 


16111.845 


4.123 


2 


642 


16115.973 


-0.326 


2 


434 


16140.839 


-8.153 


2 


730 



the mass uncertainties will generally be dominated by the 
Doppler uncertainty, focusing on low-mass stars will en- 
hance the observed signal, allowing for more meaningful 
mass constraints to be set. As proof, we again simu- 
late observations of orbiting planets, but with larger val- 
ues for K, corresponding to a less massive star or more 
massive planets. By sampling at the same times and as- 
suming the same statistical errors and jitter levels, the 
fractional error in K decreases significantly for a fixed 
number of observations. These results are also shown in 
Fig. [1] (dashed lines) . For example, a planet identical to 
Kepler-18 c orbiting a star of mass M = 0.33A/ Q would 
produce a semiamplitudc K = 15 m s™ 1 ; with only 20 
observations the RV semiamplitude could be constrained 
to within 8 percent and the stellar mass to within 30 per- 
cent. It is worth noting that these observations are all 
simulated assuming similar levels of statistical noise as 
the Kepler-18 observations. This is a reasonable approx- 
imation for the stars hosting Kepler Objects of Interest, 
but these stars are considerably fainter than the average 
Doppler planet search target. If transit timing variations 
arc detected around a considerably brighter star, as one- 
would expect from next-generation space-based planet 
finding missions, radial velocity observations could be 
carried out to considerably higher precision, decreasing 
the number of observations required to precisely measure 
the stellar radial velocity semiamplitude. 

4. SUMMARY AND DISCUSSION 

We present a method of measuring stellar and plane- 
tary masses dynamically by combining TTVs measured 
from transit light curves and follow-up radial velocity 
measurements. Our method can be used as an alter- 



0.25 



0.20 \ 

el 

S 0.15- » 



K 


= 7 m s" 1 


- - K 


=11 m s~ 


■ ■ K 


= 15 m s~ 


K 


= 19 m s~ 


k 


=23 m s~ 



Sh 

W 
13 

O 0.10 -, 

o 
td 
u 

Uh 

0.05 



°' 0< 10 20 30 40 50 60 70 80 

Total Number of RV Observations 

Fig. 1. — Derived fractional errors in the Doppler semiamplitudc 
measurement as a function of number of radial velocity observa- 
tions taken, for various values of K, the Doppler semiamplitude. 
For all observations, the same statistical uncertainties and RV jit- 
ter levels are assumed. The jitter level is 3 m s — 1 , a reasonable 
estimate for all but the youngest dwarf stars. From top to bottom, 
these curves represent semiamplitudcs of [7,11,15,19,23] m s — 1 . 
For a system like Kepler-18, where K « 7 m s , many more mea- 
surements would be required to constrain K to five percent (and 
thus the stellar mass to 15 percent). However, for a system with 
either larger planets or a smaller star, this level of precision could 
be reached with fewer observations. 



native to relying on stellar evolutionary models, which 
can be poorly constrained for non-solar type stars such 
as M-dwarfs, subgiants, and F stars. By analyzing the 
Kepler-18 system and confirming our expressions with 
dynamical simulations of this system, we show the po- 
tential of our method. 

While we show our method to be viable, especially 
for low-mass stars, using our method requires a some- 
what specific set of circumstances. The system must 
contain two planets with masses large enough to force 
a detectable Doppler signal and observable transit tim- 
ing variations on circular orbits near a first-order com- 
mensurability. Kepler data suggests planets near reso- 
nance are common: more than 12 percent of planet sys- 
tems show evidence for detectable transit timing varia- 
tions (|Ford et al.ll2012bl ). and dozens of pla nets near reso- 
nance have been confirmed through TTVs (jSteffen et al.l 
I2012bl) . Both the TTV and Doppler signals can be mea- 
sured for super-Earth planets with periods less than 30 
days; sho rt-period systems su ch as these are extremely 
common (jHoward et al.l [2012D . Thus, it is likely that 
despite the specific requirements needed to use our sys- 
tem, it can be applied to a considerable number of Ke- 
pler planetary systems. 

As shown in Equations 14 and 15 of L12, the amplitude 
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TABLE 3 

Mass Estimates for Kepler-18 system 



Object 


Cll 


L12 


Analytic Rcsult a 


Dynamical Estimate b 


Star (M ) 
Planet c (M e ) 
Planet d (M e ) 


0.972 ± 0.042 

17.3 ± 1.8 

16.4 ± 1.4 


Assumed Cll 
20.2 ± 1.9 
17.4 ± 1.2 


0.83 ±0.51 
18.6 ± 11.6 
15.4 ±9.5 


Q2+°' ei 
u,y -0.40 

14 8+ 9 ' 4 

15 4+11-0 
10 -7.0 



a Result derived by applying Equations 11-13 
b Result determined from numerical integrations 



of the TTV signal is given such that 

|V|~|V dlunp ed|(l + ^), (15) 

with | Vdampedl the amplitude of the TTV signal if the sys- 
tem were damped of its free eccentricity. The quantity 
Z[ ree is defined such that Zf ICC = /zf roc + gz' ilcc , where z 
is the complex eccentricity of the planet, z = eexp(icj). 
Thus, our method as described will break down unless 
l-Zfreel <C |A|. This is a reasonable assumption for plan- 
ets with periods under ten days. In theory, even if a 
non-negligible amount of free eccentricity remains in the 
system, a detailed radial velocity orbital solution and 
light curve analysis could be used to calculate Zf ICC and 
determine the stellar and planetary masses. 

As a projection of the utility of this method, consider 
KOI 1241, a system containing two planets with periods 
of 10.5 and 21.4 day s orbiting a giant star (R = 3.14i? Q , 
ISteffen et al.l l2012bl) . From the transit timing signal, 
there is evidence that this system has not dissipated its 
free eccentricity, meaning it is not optimal for our study. 
However, with nine quarters of public Kepler data, wc 
can constrain the TTV signal caused by the larger planet 
in this system to 8.2 percent. Moreover, with only nine 
radial velocity observations, the radial velocity semi- 
amplitude of the larger planet can be measured to 5.3 
percent. Thus, from our method alone, if a system ex- 
isted that was nearly identical to KOI 1241 but damped 
of free eccentricity, by Equations 11-13, we expect we 
could determine the stellar mass to 29 percent and the 
mass of the larger planet to within 22 percent. Our 
method could also be applied to KOI 1241 in the future 
if enough radial velocity data is collected to determine 
the magnitude of Zf ICC for the system. The uncertain- 
ties in the TTV signal of KOI 1241 are larger than the 
uncertainty in the RV semiamplitude. The transit tim- 
ing errors will decrease as more Kepler data are released: 
decreasing the TTV error to five percent without includ- 
ing any additional radial velocity observations will reduce 
the uncertainty in the stellar mass to 20 percent. Thus 
our method could provide significant constraints on stel- 
lar masses in regimes where stellar atmospheres are less 
well-understood, such as subgiants and cool stars. For 
these cases, our method will be able to compliment as- 
teroseismology results as an independent measure on the 
mass of the star. Moreover, our method can be used 
to find systems where the analytic stellar mass is sub- 
stantially different than the Kepler Input Catalog values, 
which can then be followed up with dynamical modeling, 
asteroseismology, or high resolution spectroscopy to bet- 
ter characterize the star and orbiting planets. 

Our method of estimating the stellar mass of a planet 



host star requires a many high-precision radial velocity 
measurements as well as observations of enough measure- 
ments to measure TTVs. Traditionally, for radial veloc- 
ity searches the stellar mass is estimated by interpolat- 
ing spectroscopic properties such as the stellar effective 
temperatu re and metallicity on to stellar evolution model 
grids (e.g. iSozzetti et al.ll2007l ). From a careful analysis 
of transit light curves, one can measure the ratio a/R* 
and the duration of the transit event, which provide an 
estimate of the stellar density. However, stellar evolution 
models arc still required to break the degeneracy between 
stellar radius and mass. While these traditional meth- 
ods allow for precise measurements of stellar masses, the 
results are mo del- dependent and may be subject to sys- 
tematic errors induced by the stellar models. This is 
particularly true for stars very different from the Sun. 
Our method has the advantage of providing an accurate, 
model-independent stellar mass, albeit with a trade-off 
in precision. However, even in cases for which the preci- 
sion of our method is poor, the resulting mass estimate 
can be used as a prior in the more traditional model-grid 
interpolation. 

With present data our technique is only viable as an 
alternative to stellar modeling in the most exceptional 
cases. Transit timing variations have been detected to 
remarkable precision by Kepler, but very few KOIs have 
been followed up with radial velocity measurements. In 
the cases where RV data exists, only enough measure- 
ments were collected to confirm the planetary nature of 
the syst em, not to independ ently measure the planetary 
masses (jHolman et al.1 [20101 ) . Our routine will become 
more useful for systems around brighter stars, when more 
RV measurements can be efficiently collected and the RV 
semiamplitude can be better constrained. 

Despite the faintness of the Kepler planet candidate 
host stars, there are a few stars that would be ideal can- 
didates for applications of our method. From the collec- 
tion of Kepler Objects of Interest, we searched for stars 
hosting at least two transiting planets each with P < 25 
days. We required at least one planet to be larger than 
2 i?^ and the planet periods to lie within five percent 
of a first-order mean-motion resonance. To ensure that 
all targets were optimized for radial velocity follow-up, 
we eliminated all targets fainter than mx p = 13.0. After 
making these cuts, wc find 8 candidate systems to which 
this technique can be applied: KOI 85, 111, 115, 117, 
244, 304, 1241 and 1930. As stated earlier in this section, 
KOI 1241 is not an ideal target because it has not been 
fully damped of its primordial eccentricity and there is 
not enough radial velocity information to uniquely deter- 
mine the eccentricity of both planets. Of the remaining 7 
systems, the CPS team has collected more than 10 radial 
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velocity measurements only on one, KOI 244. Additional 
radial velocity measurements of any or all of the above 
systems would enable further validation of our procedure 
as well as additional constraints on the masses of each of 
the stars and their planets. Moreove r, next-generation 
plane t finding missi ons, such as TESS rown fc Latham! 
[20081) and PLATO (jCatala et al.|[2010t) will target bright 
stars, making detailed radial velocity follow-up observa- 
tions of systems exhibiting transit timing variations a 
much more practical possibility. 
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TABLE 4 

Transit Times for Kepler Transiting Planet Candidates 



KOI 


n 




TTV„ 








BJD-2454900 


(d) 


(d) 


137.01 





198.3142 


-0.0006 


0.0012 


137.01 


1 


205.9557 


0.0002 


0.0013 


137.01 


2 


213.5973 


0.0019 


0.0018 


137.01 


3 


221.2389 


0.0017 


0.0010 


137.01 


i 


228.8804 


0.0021 


0.0017 


137.01 


5 


236.5220 


0.0014 


0.0015 


137.01 


6 


244.1636 


0.0031 


0.0011 


137.01 


7 


251.8052 


0.0030 


0.0011 


137.01 


8 


259.4467 


0.0037 


0.0012 


137.01 


9 


267.0883 


0.0050 


0.0012 


137.01 


10 


274.7299 


0.0041 


0.0012 


137.01 


11 


282.3714 


0.0037 


0.0018 


137.01 


12 


290.0130 


0.0035 


0.0018 


137.01 


13 


297.6546 


0.0028 


0.0012 


137.01 


14 


305.2961 


0.0034 


0.0011 


137.01 


15 


312.9377 


0.0026 


0.0011 


137.01 


16 


320.5793 


0.0013 


0.0011 


137.01 


17 


328.2209 


0.0018 


0.0012 


137.01 


18 


335.8624 


0.0008 


0.0015 


137.01 


20 


351.1456 


0.0001 


0.0011 


137.01 


21 


358.7871 


-0.0005 


0.0012 


137.01 


22 


366.4287 


-0.0027 


0.0014 


137.01 


23 


374.0703 


-0.0033 


0.0014 


137.01 


24 


381.7119 


-0.0024 


0.0011 


137.01 


25 


389.3534 


-0.0025 


0.0011 



8 



Montet and Johnson 

TABLE 4 — Continued 



KOI 


n 




TTV„ 




& n 








BJD-2454900 


(d) 




(d) 


137. 


01 


9f; 


396 


9950 


n 00*3*3 


u 


001 1 
UU11 


137. 


01 


27 


404 


6366 


n oo*3fi 





0013 


137. 


01 


28 


412 


2781 


-U.UU1U 


o 


0012 


137. 


01 


20 


419 


.9197 


-0.0037 


o 


0011 


137. 


01 


30 


427 


5613 


n oo*38 

-u.uuoo 


o 


0015 


137. 


01 


31 


435 


.2028 


o oo*3Q 

-u.uuoy 


o 


0013 


137. 


01 


32 


442 


.8444 


-0.0019 


o 


0012 


137. 


01 


33 


450 


4860 


-0.0013 


o 


0011 


137. 


01 


34 


458 


1276 


-0.0017 


o 


0010 


137. 


01 


35 


465 


7691 


-0 0008 

u. uuuu 


o 


0015 


137. 


01 


36 


473 


4107 


n ooo^ 

U . UUUU 





0013 


137. 


01 


37 


481 


0523 


0.0016 


o 


.0011 


137. 


01 


38 


488 


6938 


U . UUiU 





0014 


137. 


01 


39 


496 


3354 


0.0027 


(") 


0015 


137. 


01 


10 


503 


9770 


o 00*30 

u . uuou 





.0014 


137. 


01 


44 


511 


6186 


00*38 
u . uuoo 


("J 


0011 


137. 


01 


42 


519 


2601 


00*39 
u . uuo^ 


(") 


0013 


137. 


01 


43 


526 


9017 


00^9 

U . UUJ^ 





.0013 


137. 


01 


44 


534 


5433 


n nn^ 
u . uuoo 


("J 


0013 


137. 


01 


45 


542 


.1848 


0.0034 





.0010 


137. 


01 


16 


549 


8264 


009*^ 
u . uu^o 





.0014 


137. 


01 


47 


557 


4680 


004Q 
u . uuiy 


("J 


0021 


137. 


01 


18 


565 


1095 


00*3^ 
u . uuoo 





.0011 


137. 


01 


49 


572 


7511 


00*30 
u . uuou 





.0012 


137. 


01 


50 


580 


3927 


0.0017 


("J 


0010 


137. 


01 


51 


588 


.0343 


0.0034 





.0016 


137. 


01 


52 


595 


.6758 


0.0013 





.0011 


137. 


01 


54 


610 


9590 


n ooor 

-u . uuuu 





0012 


137. 


01 


55 


618 


6005 


-0.0019 





0012 


137. 


01 


56 


626 


2421 


-0.0015 





0016 


137. 


01 


57 


633 


8837 


n nn9fi 

-u.uu^u 





0012 


137. 


01 


58 


641 


.5253 


-0.0016 





0012 


137. 


01 


61 


664 


4500 


-u.uuiu 





0014 


137. 


01 


(] ■? 


672 


,0915 


-U . UUiiO 





0012 


137. 


01 


63 


679 


.7331 


-0.0042 





0012 


137. 


01 


65 


695 


0162 


00*39 
-u.uuoo 





0011 


137. 


01 


66 


702 


,6578 


-0 009^ 





.0013 


137. 


01 


67 


710 


2994 


-0.0029 





0012 


137. 


01 


68 


717 


9410 


009Q 
-u.uu^y 





0011 


137. 


01 


60 


725 


.5825 


000/L 

U. UUU'i 





0013 


137. 


01 


71 


740 


8657 


OOOfi 
u . uuuu 





.0013 


137. 


01 


72 


748 


5072 


OOOfi 
u . uuuu 


(") 


0014 


137. 


01 


73 


756 


1488 


0.0014 





.0011 


137. 


01 


7-1 


763 


.7904 


004Q 
u . uuiy 





.0017 


137. 


01 


75 


771 


4320 


0.0011 


("J 


0016 


137. 


01 


76 


779 


0735 


00*3*3 
u . uuoo 





.0012 


137. 


01 


77 


786 


7151 


0.0021 


("J 


0011 


137. 


01 


78 


794 


3567 


0.0031 


(") 


0012 


137. 


01 


70 


801 


9982 


00*3Q 
u . uuoy 





.0012 


137. 


01 


80 


809 


6398 


009Q 
u . uu^y 


("J 


0014 


137. 


01 


81 


817 


2814 


00*38 
u . uuoo 


(") 


0011 


137. 


01 


89 


824 


9229 


0.0016 





.0012 


137. 


02 


Q 


194 


,8832 


0000 
u . uuuu 


(") 


0010 


137. 


02 




209 


,7421 


-0.0012 


(") 


0011 


137. 


02 


2 


224 


,6011 


009"^ 





0010 


137. 


02 




239 


,4600 


-0.0037 





0010 


137. 


02 


4 


254 


,3189 


-0.0029 





0010 


137. 


02 


ji 


284 


,0368 


00*39 
u.uuoz 





0011 


137. 


02 


7 


298 


,8958 


0040 

-U .UU4U 





0011 


137. 


02 


g 


313 


,7547 


-0.0013 





0010 


137. 


02 


9 


328 


,6136 


009*3 





0023 


137. 


02 


10 


343 


4726 


000 ^ 
-u. uuuo 





0011 


137. 


02 




358 


,3315 


0008 
u . uuuo 





.0009 


137. 


02 


12 


373 


1905 


0.0017 





.0012 


137. 


02 


13 


388 


,0494 


00*3fi 
u . uuou 


("J 


0010 


137. 


02 


14 


402 


.9083 


0.0029 





.0011 


137. 


02 


15 


417 


,7673 


0.0020 





.0010 


137. 


02 


16 


432 


,6262 


0.0008 





.0012 


137. 


02 


17 


447 


,4852 


0.0008 





.0011 


137. 


02 


18 


462 


,3441 


0.0008 





.0011 


137. 


02 


19 


477 


,2030 


-0.0001 





.0011 


137. 


02 


20 


492 


.0620 


0.0002 





.0010 


137. 


02 


21 


506 


,9209 


-0.0023 





.0011 


137. 


02 


22 


521 


,7799 


-0.0024 





.0010 


137. 


02 


23 


536 


.6388 


-0.0038 





.0010 
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KOI 


n 




TTV„ 




a n 






BJD-2454900 


(d) 




(d) 




21 


kc:i /in 77 


n oo*^ 


n 
u 


001 


137.02 


25 






o 


0011 


137.02 


26 


^Kl 91 ^fi 


n oo 1 ^ 


o 


0010 


137.02 


27 




-0.0013 


o 


0015 


137.02 


29 


fi9^ 7Q94. 




o 


0010 


137.02 


31 


fi 1 ^ ^1 03 


0039 


Q 


.0014 


137.02 


32 


U t U.OUSO 


0.0015 


Q 


.0011 


137.02 


33 


fi&^ 9989 


009fi 


g 


0009 


137.02 


31 


700.0871 


-0.0011 





.0011 


137.02 


35 


714.9461 


0.0008 





.0010 


137.02 


36 


729.8050 


-0.0019 





.0010 


137.02 


37 


744.6640 


-0.0012 





.0011 


137.02 


38 


759.5229 


-0.0027 





.0011 


137.02 


39 


774.3818 


-0.0030 





.0012 


137.02 


40 


789.2408 


-0.0034 





.0010 


137.02 


11 


804.0997 


-0.0052 





.0010 


137.02 


12 


818.9587 


0.0003 





.0013 



